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Entanglement in incoherent mixtures of pure states of two qubits is considered via 
the concurrence measure. A set of pure states is optimal if the concurrence for any 
mixture of them is the weighted sum of the concurrences of the generating states. 
When two or three pure real states are mixed it is shown that 28.5% and 5.12% of 
the cases respectively, are optimal. Conditions that are obeyed by the pure states 
generating such optimally entangled mixtures are derived. For four or more pure 
states it is shown that there are no such sets of real states. The implications of these 
on superposition of two or more dimerized states is discussed. A corollary of these 
results also show in how many cases rebit concurrence can be the same as that of 
qubit concurrence. 
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Ctf ■ I. INTRODUCTION 

cr 

Entanglement properties of pure and mixed quantum states have been the subject of in- 
^ \ tense and extensive study in the recent past Of these, entanglement in qubits or spin- 1/2 
systems have dominated due to their use as fundamental objects in quantum computations. 
For an arbitrary state of two qubits the concurrence measure (or its square, called tangle) 
introduced by Hill and Wootters ^ is simply calculable from the density matrix and is a 
measure of entanglement. To be precise the entanglement of formation [3] is a monotonic 
function of the concurrence. The concurrence measure has been extensively applied in many 
physical contexts, for instance in the study of quantum phase transitions [4J . In a collection 
of qubits, concurrence measures the entanglement present within any chosen pair. Thus due 
to the monogamy property of entanglement [B[ it is reasonable to expect that states with 
large multipartite entanglement have low or vanishing concurrence. In fact for a random 
state with more than six qubits the probability that a chosen pair has nonzero concurrence 
is vanishingly small, most of the entanglement is of the multipartite kind 0, 0] ■ 
To elaborate on this property, consider a mixed state of two qubits 

k 

p = hkXV'il* k>i, (i) 

where pi > 0, Y^iVi — 1 an d the projectors are arbitrary, in particular, \ipi){ipi\ need not be 
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orthogonal. The convexity of concurrence |8j implies that 

k 

C(p) < X>C7(iM (2) 



i=i 



where C is the concurrence function, an entanglement monotone Thus the maximum 
that C(p) can attain is the weighted sum of the concurrence of the extremal (pure) states. If 
there exists a set (IV'i)) ■ • ' j^fe)} sucri that equality is obtained in Eq. (j2J) for any arbitrary 
set of weights {pi}, it is referred to herein as optimal. However note that such a property 
will be specific to the concurrence measure of entanglement. 

For states that are real in the standard basis, it is shown that a very large fraction 
of states made by incoherently superposing 2 two-qubit states optimize their entanglement. 
This property is analyzed in detail in this paper and conditions to be satisfied by the extremal 
states such that the resultant density matrix is optimal are derived. Any real density matrix 
can be tested for optimality of its diagonal decomposition using the inequalities derived. 
These are also generalized beyond two states, and it is shown that for more than three real 
states, not one optimal decomposition exists. 

The relevance to superposed dimers will be studied in section IHIl The relation between 
entanglement and superpostion of quantum states is an interesting one [||. Indeed super- 
position of states with a tensor product structure is necessary for entanglement, however of 
course this is not sufficient. There is a significant amount of literature establishing bounds 
on various entanglement measures for the superposition in terms of the entanglement in the 
states that are being so superposed 0-15 . 



II. RANK-WISE STUDY OF OPTIMIZING MIXTURES 

If two pure and real states and l^) are chosen at random, it is shown in this 
paper that in 28.5% of cases the resultant entanglement in p\ipi)(ipi\ + (1 — p^V^XV^I is 
the maximum possible, namely equality holds in Eq. (J2J). This implies that on superposing 
2 two-qubit states, 28.5% of the states will remain entangled optimally, as defined in the 
introduction. The fraction 0.285 of optimal pairs is interesting and strong evidence that it 
is actually (it — 2)/4 is presented in an Appendix. 

The conditions under which such optimality occurs is now obtained, however this is done 
in a more general setting in what follows. In particular, extending these results to arbitrary 
mixtures of three real pure states, one finds that in about 5.1% of cases this gives rise to 
optimally entangled states. It is also then shown that for four or more states there is not even 
one set of pure real states, such that all their mixtures are optimal. These generalizations 
are of relevance when more than 2 two-qubit states are superposed. 

The reader is first reminded of the procedure to find the concurrence in p, a given state 
of two qubits 0. The spin-flipped state p = a y ® o y p*o y ® a y is found, where the complex 
conjugation is done in the standard basis. Then the matrix pp is diagonalized and has 
positive eigenvalues \i\ > /i 2 > /x 3 > The concurrence C(p) is max(0, yZ/Tj" — y 7 /^ — 

This somewhat involved definition of the concurrence renders it opaque for considerations 
of optimality. However it is possible to express the concurrence of p = ^*L 1 Pi|V , «) (V'il more 
explicitly in terms of the the states and the weights Restrict to the case k < 4, that 
is the size of the generating set of states is not larger than the maximum rank of p. Note 
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that if any p is expressed in its eigenbasis this is not a restriction at all. It is now shown 
that the eigenvalues of pp are the same as that of r'r'* where 

r ij = VPiPj r *J> aIld r V = (i>i\ a V ® <T y |^>. (3) 

Thus rather than using the density matrix directly, the pure states comprising a particular 
ensemble are used. Note that |rn| and 1 7^22 1 are the concurrences, C(V>i) and C^fe), of the 
pure states \ipi) and \ip 2 ) respectively. 

For convenience the eigenvalue equation for pp is considered, whose right eigenvectors 
can be written in the nonorthogonal, sub- normalized basis of the extremal states as \u) = 
Ylii a iWi)i where |^-) = y/pi\ipi}. Also, writing p = J2i=i 1^0(^1 and using the fact that 
(i>'m\PPW) = nfy'mW) results in 

k k 

P'm^'ji ~ PKni)<Xi = 0, l<m<k, (4) 

i=l j=l 

where 

t'ji = y/PjPi tji = y/PjPi(lpj\lpi) (5) 

is a matrix of inner products (the Gram matrix) and 

P'mj = VP^Pj Pmj = ^P m Pj{^m\p\^j}- (6) 

Note also that since, p* = J2Pi\i J i)(i'i\i 

P' mj = $m\pWj) = (tfmWv ® VyC^PiWi) {^i\)°y ® <r y Wj) (7) 

the following matrix identity is readily derived: p' = r'r'* = r'r'\ where r and r' are defined 
above in Eq. (J3J). 

For the k equations in Eq.(jH) to have non-trivial solutions det(A) = 0, where 

A m j = (p't')mi — p(t') m i, (8) 

which further implies that 

det(t') det(rV* -/i/) = 0. (9) 

As the number of vectors \ipi) are no larger in number than the dimensionality of the Hilbert 
space, we assume them to be independent and therefore det(t') = p\ - ■ -pk det(t) 7^ 0. Thus 

det(rV* - pi) = 0. (10) 

and hence the characteristic polynomials of pp and r'r'* are identical. 

If the state p is real in the computational basis then the eigenvalue problem of r'r'* is 
that of r' 2 , whose eigenvalues are the square of the eigenvalues of r', which are indicated 
as A. The expression for concurrence is derived and the conditions for optimality are now 
considered case-by-case starting from k = 2. 
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A. Rank-2 density matrices: k = 2 

For the case of mixtures of two real pure states, k = 2, the above considerations lead to 
the following characteristic equation of r' (as defined in Eq. OH])): 

A 2 -£A+2Wi2 = (11) 

where £ = p±r u + p 2 r 22 , and X12 = ( r n r 22 — r i 2 )- The eigenvalues are therefore 

A± = ^(e± ^-^2X12) (12) 

If X12 > and f > then A + > A_ > 0, and C(p) = A + - A_. Alternatively if X12 > 
and f < then A_ < A + < 0, and C(p) = |A_| — |A + |. Thus if X12 > 0, irrespective of the 
sign of £, we have that 

c( P ) = ve-4 PlP 2Xi2 < iei < PiCtyi) + P 2C(^ 2 ), (13) 

confirming the convexity of concurrence. 

The case X12 < is more interesting and leads to optimality. If £ > then A_ < < A + 
and A + > |A_|. It follows that C(p) = A + — |A_| = A + + A_ = f . Combining a similar 
analysis of the case f < one gets that when X12 < 0, irrespective of the sign of f 

C(p) = |f I = bir n +p 2 r 22 |. (14) 

This brings us to the possibility that if r u r 2 2 > then C(p) = PiC^^+^C^)- Thus when 
l^i) and 1^2) satisfy the conditions that r u r 2 2 > and rur^ — r\ 2 < any arbitrary mixture 
of these pure states has the maximum entanglement which is their average entanglement. 
This is the first set of optimality conditions that we derive. 

The set of such optimal states is a subset from pairs of real states. Each real state of two 
qubits is characterized by 4 real coefficients, say X{, i — 1, . . . , 4. The normalization condition 
means that there is a isomorphism between these and the 3-sphere x\ + x\ + x\ + x\ = 1. 
Apart from the fact that states differing by a sign are really the same (thus the space is a 
projective space) the states maybe thought of as points in S* 3 . Thus a pair of real states is 
a point on the manifold S 3 x S 3 , and the set of optimal states forms a subset therein whose 
fractional volume is of natural interest. 

Assume that the real states of two qubits are distributed uniformly on S 3 , namely choose 
the Haar measure. Equivalently, the probability density of random real pure states [l6| of 
two qubits is given by 

where Xt are the state components in a generic basis such as the computational one. The 
fraction of optimal states f 2 is then the following integral on S 3 x S* 3 , written in terms of 
the ambient space components (xj, yi) of 1R 4 x R 4 : 

h = I ©(rnr 22 )e(r? 2 - r u r 22 )P({x t })P({ yi })l[dx t dy t . (16) 
•* i=i 
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Here r n = (V>i|<7» <g> <r v |^) = 2(x 2 x 3 - Xi£ 4 ), r 22 = (V'aky ® 0^2) = 2(^2/3 - 2/12/4), 
r i2 = (V'i I ® = ^22/3 + 2/2^3 — ^i2/4 — ^iUi, an d © is the Heaviside step function that 

is 1 if its argument is positive and zero otherwise. An exact evaluation of this integral seems 
possible and equal to (tt — 2)/4 (see the Appendix), however it is quite easy to simulate the 
process by choosing two independent vectors Xi and yi distributed according to Eq. (TT5l) and 
checking to see if the optimality condition is satisfied. The initial choice of vectors is done by 
simply taking 4 numbers from any zero centered normally distributed set and normalizing 
them. This procedure gives the fraction / 2 to be approximately 0.285, in good agreement 
with the value of (it — 2)/4. 



B. Rank-3 density matrices: k = 3 

Now we consider the general setting of mixing three real pure states, k = 3, which leads 
to the cubic equation for the eigenvalues of r', (defined in Eq. Q): 

Z(A) = A 3 + £iA 2 + 6A + £ 3 = 0, (17) 

where the coefficients £j are 

6 = - ^Piru, £2 = ^PiPj^iiTjj - rj), £ 3 = -P1P2P3 det(r). (18) 

i ijtj 

Now we state two Lemmas which are key to understanding the nature of the roots of cubic 
equations and the possibility of optimal states in the case k = 3. 

Lemma 1. If p(x) is a cubic in x, p(x) = x 3 + ax 2 + bx + c with real coefficients, has real 
roots, and is such that a, b, c < then p(x) = has two negative roots and one positive root 
with the positive root being greater than the other two in modulus. 

Proof. Since c < the product of the roots is positive which implies that either all roots are 
positive or there are two negative and one positive root. However since b < the quadratic 
p'(x) has one positive and one negative root. Hence all the roots of p(x) cannot be positive 
since the roots of p'(x) have to lie between the roots of p(x) by Cauchy's mean value theorem 
for differentiable functions. Hence the polynomial has two negative and one positive root. 
Observe that a < implies the sum of the roots is positive which implies that the positive 
root has the largest modulus. □ 

Lemma 2. If p(x) be a cubic in x, p(x) = x 3 + ax 2 + bx + c with real coefficients and has 
real roots, and is such that b < and a, c > then the cubic has two positive roots and one 
negative root with the negative root being greater than the other two in modulus. 

Proof. The proof of this lemma is on similar lines to the previous one. □ 

Case 1. Now suppose we have £1, £3 < then along with the condition that £ 2 < we have 
by Lemma 1 that the characteristic equation has two negative and one positive root and that 
the positive root has the largest magnitude. Without loss of generality let us assume that 
Ai > then we have the square-roots of the eigenvalues of pp = (a y ®o- y p) 2 are Ai, — A 2 , —A3. 
Hence the concurrence of p is C(p) = max(0, Ai + A 2 + A3) = max(0, — £1) = |£i|. Thus 



C(p) = pir u + P2T22 + P3^33- 



(19) 
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If we have > for i = 1, 2, 3 then clearly C(ipi) — r %%- Also if £ 2 < then since p^s can be 
arbitrary positive reals bounded by 1 we have each term PiPj{rnTjj — r? ) < which means 
ruTjj — rfj < 0, i 7^ j is a necessary condition. Thus we have that if 

r u > 0, r u rjj - rj < 0, i ^ j (20) 

and 

det(r) = ^Vn(r 22 r 33 - r\ z ) - 2(rnr 22 r 33 - ri 2 r 23 r 3 i) > (21) 

eye 

then any mixture of the triple {|?/>i), {ipz), 1^3)} will be optimally entangled, that is 

3 

C(p) = 5>C7(ifc) (22) 

Case Similarly suppose we have £1,^3 > then along with the condition that £ 2 < 
we have by Lemma 2 that the characteristic equation has two positive and one negative root 
and that the negative root has the largest magnitude. Let us assume that Ai < then we 
have that the square-roots of the eigenvalues of pp = (a y <S> cr y p) 2 are — Ai, A 2 , A3. Hence the 
concurrence of p C(p) = max(0, — (Ai + A 2 + A3)) = max(0,£i) = £1, since £1 > 0. Thus we 
have if 

ra < 0, ruTjj - r?- < 0, i ^ j (23) 

and 

det(r) = ^ r n( r 22^33 - rl 3 ) - 2(riir 22 r 33 - ri 2 r 23 r 31 ) < (24) 

eye 

then again any mixture of the triple {IV'i)? IV^), \i J 3)} is optimally entangled. These condi- 
tions that are derived for optimality of entanglement for all mixtures are both necessary and 
sufficient. They can be compactly stated as, 

r«det(r) > 0, [r]« < 0, i = 1,2,3, (25) 

where [r]a are the principal minors of r. 

Once again we can therefore write the fraction of triples of pure real states of two qubits 
whose arbitrary mixtures are optimally entangled in terms of an integral such as in Eq. (1161) . 
which when evaluated using a procedure extending the previously described one, gives the 
fraction ^3 = 0.0512 or 5.12%. Parenthetically, while it is possible that this integral in all 
its complexity evaluates exactly as well, the authors are fairly confident that they cannot 
find it. It is interesting that indeed optimal triples are such that any two of them are also 
optimal. This has to be the case as the vanishing of any one of the probabilities pi must 
still be optimal. Thus when we go from mixtures of rank-2 to rank-3 we see a drastic drop 
in the percentage of states that lead to optimal entanglement. 

C. Rank-4 real density matrices are never optimal: k > 4 

Incoherently superposing four of more pure and real states leads to a qualitatively different 
behavior, as shown below. When k = 4 the rank of the eigenvalue problem for pp is full, 
in the sense that it is the dimensionality of the Hilbert space. From our discussion above 
it is clear that we have a quartic polynomial whose constant term is det(r). For optimality 
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we have to have the concurrence evaluating to the trace of r', we need to have either one 
eigenvalue positive and three negative eigenvalues, all of them smaller than the positive one 
in modulus, or one negative eigenvalue higher in modulus than three positive eigenvalues. 
In either case this implies that det(r) < 0. However it is easy to prove that det(r) > 0. 
Indeed, 

r\j = (ipi\cr y <8> o- y \ipj) = 22 Fim( m \°~y ® v y \n)Fnj- (26) 



mn 



Here F im = (^|m), and \m) is any real orthogonal basis, for instance it could be the 
computational one. Therefore 

det(r) = det[a y <g> a y ] det(F) 2 = det(F) 2 > 0, (27) 

the final inequality following from the reality of the transformation functions. Therefore 
unlike the rank-deficient cases the sign of det(r) is always positive and this rules out the 
existence of even one real quadruplet such that any arbitrary mixture of these remains 
optimally entangled. This obviously implies the non-existence of even one set of real optimal 
state for k > 4. It is necessary to have complex states in the ensemble for optimizing the 
entanglement in this case. 

Entanglement in real qubits have been studied earlier by also restricting the Hilbert 



space to the space of reals, the so-called case of "rebits" [17J. In this case minimization 
of the entanglement is also carried out only over the real ensembles that are realizations 
of the density matrix, unlike in the current paper, where we have used the usual formula 
for concurrence. The rebit formula for concurrence is \tr(a y <g> <J y p)\ which in terms of the 
quantities introduced in this paper is | YliiPi r a\- ^ n the case when k = 2 this is also the 
actual concurrence if we only require that X12 = ^11^22 — i\i < 0, relaxing the condition 
that r u r 2 2 > 0. This is true in about 78.5% (7r/4 fraction) of cases when k = 2. Thus 
stated in terms of rebits, the present work also implies that for 78.5% of pairs of real pure 
states that are mixed, the rebit entanglement coincides with the usual qubit entanglement. 
Similar generalizations to the case of triples of states gives about 48.5%, which results from 
relaxing the conditions that all Ta have the same sign. However it follows from the above 
considerations that when we take a mixture of four (or more) real states, the resultant rebit 
entanglement is always suboptimal to the usual concurrence obtained on the full complex 
Hilbert space. 



III. DIMERIZED STATES AND OPTIMALITY 

As an application of the study of optimizing mixtures of two qubits, the problem of 
entanglement sharing in superpositions of states with a dimerized structure is now taken 
up. If there are many pure states of 2N qubits such that qubits 1 and 2 are entangled only 
with each other, 3 and 4 with each other and so on, and each of these pairs are in pure 
states, then an implication of the previous section is that superposing such "dimerized" 
states results in rather robust entanglement, especially if only two or three such states are 
superposed. On adding more such states, the entanglement in the pairs comes down due to 
the lack of optimality, and will lead to more global, or multipartite entanglement. Consider 
the, in general unnormalized, state: 

k 

lV>>=X>®£il#>; where m=^\00) J +^\01) J+1 }\10) J + 5}\ll) J , (28) 

i=l 
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and ^ | a^l 2 = 1. Thus the state is a superposition of k states, labelled by i, each of 
which has N pairs of entangled two qubit (normalized) pure states, labelled by j. No 
two pairs are entangled with each other. Such superpositions arise in many context, for 
example in the Resonating Valence Bond states [l8|. However note here that the "dimers" 
superposed are of the same kind, that is, the entangled pairs of particles are the same. Most 
Hamiltonian systems have some form of time-reversal ( ant i- unitary) symmetry that renders 
their eigenstates real. 

Henceforth the state of the first entangled pair in the i th state is referred to as 
For simplicity consider the case when k = 2. Let a x = cos(#) and a 2 = sin(#). As the 
superposed states are not orthogonal, there is the normalization factor A/", for the state 
in Eq. (128]) which is 



N 

N = 1/^/1 + sin 20//i, where m = \{{a)a] + + 7)7? + 5)5]). (29) 

The reduced density matrix of any two qubits that are entangled in the original states 
which without loss of generality can be taken as the first two qubits, is 

p 12 =M 2 (cos 2 0|^i)(^i| + sin 2 ^|^ 2 )(^ 2 | +/i 2 sin0cos0(|^i)(^ 2 | + • (30) 
Here p 2 is defined as 

N 

iii = IlKH 2 + + iH + 5 ) 5 D- ( 31 ) 

i=2 

In all generality this is all that can be said about P12, however for most states the interfer- 
ence or coherence term is negligible, due to the typical smallness of /ii and /12. Thus the 
approximation p i2 ~ p, where 

p = cos 2 ^)^! + sin 2 0|^ 2 )(^ 2 |, (32) 

is a good one. 

To estimate the typical value of p\ and p 2 consider the Hilbert space of each entangled 
pair consisting of two qubits and take for the distribution of the coefficients the uniform 
or Haar measure of Eq. (fT5l) . The averages of p\ and p 2 are both zero in this ensemble. 
However the second moments are nonzero and can be shown to be 

This follows on observing that p\ and p 2 in Eqs. (|29| [3TT) are N and N — 1 fold products of 
inner-products of two four- vectors. Consider one such (square of the) inner-product and its 
ensemble average, the average being over the distribution where each of the 4-vectors, here 
denoted simply as Xi and i/i are distributed according to the measure in Eq. (TT5]) : 




= 1/4. (33) 



i=l 



The first equality follows as the odd powers average to zero, and the second equality follows 
as the ensemble average of each of the xj is 1/4, which follows most easily from normalization. 
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FIG. 1: Concurrence of the first two qubits vs p\ = cos 2 (6 l ) for the superposition of two 
states {k = 2) with N = 4 and N = 8 (see Eq. fl2"B"l) ) for two different realizations (left and 
right columns). The incoherent case is obtained on neglecting the interference terms in the 
superposition as in Eq. (132]) . The bottom left panel shows a case of optimality. 



Thus /ii and \i2 are typically of the order of 2~ N and 2~^ N ~^ respectively. In practice 
for most states with N > 8 it is difficult to distinguish whether the incoherent two-qubit 
state p obtained on dropping the interference term as in Eq. (132]) is used or the actual 
reduced density matrix p 12 (Eq. (13T)]) ) is used. This is illustrated in Fig. (pQ), where random 
realizations of the two-qubit states |0* ) are used in Eq. (1281) with k = 2. Two such realizations 
are selected for the cases of N = 4 and N = 8 and the concurrence in the reduced density 
matrix of the first two qubits are calculated based on the exact state pu (Eq. fl30|) . referred 
to in the figure as "coherent") and the approximation p (Eq. (1321) . referred to in the figure 
as "incoherent"). The concurrences are plotted as a function of the mixing between the 
states that are superposed in Eq. (128]) . It is seen from the figure that when N = 8 whatever 
difference persist between the entanglement in p and p^ is not visible, while it is for N = 4. 
Also from the N = 8 case, the realization on the right, illustrates the convexity of concurrence 
as at p = 0, 1 two different two-qubit pure states are obtained while for intermediate p, the 
density matrix is an incoherent superposition of these (see Eq. (132]) ). However for N = 8 
the realization on the left is peculiar in that the concurrence for the intermediate values of 
p is just a linear interpolation of the concurrences of the pure states. In this case one has 
hit upon what is studied in this paper as an optimal pair of pure states. 

The approximate form of the reduced density matrix in Eq. (1321) obtained by using the 
incoherent superposition condition holds also for complex states, however optimality results 
arise in the case of real state as studied in detail in the previous section. For k > 2 similar 
considerations give rise to states of the form as in Eq. ([T]) where pi are equal to |cij| 2 . In terms 
of dimerized states, we see that superposing two of them leads to pairs of qubits that were 
originally entangled with each other retaining much of it. About 28.5% of pairs of qubits 
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would have simply the weighted entanglements of the states before superposition. On the 
other hand superposing three dimerized states leads to a significant decrease in the fraction 
of robust dimers, which is now 5.12%, and superposing more than three, the entanglement 
of a pair is bound to be smaller than the weighted entanglements. 



This work has used a particular measure of entanglement of two-qubit density matrices, 
namely concurrence, and studied questions of optimality, as defined herein, mainly restricting 
to the space of real states. Concurrence is only one measure of entanglement, but the 
entanglement of formation being a monotonic function of it renders it rather unique and as 
such this measure has been used in very many studies. What exactly the issue of optimality 
says about the geometry of the quantum space of states [19f is an interesting question that 
is not pursued here. It may also be interesting to study other measures of non-classical 
correlation, such as discord [20| from this perspective. 



One possible application of optimality was discussed in the superposition of dimerized 
states. Such states are found in many quantum spin systems, such as the Majumdar Ghosh 
hamiltonian j2l[, and superpositions are relevant in the neighborhood of avoided crossings 
where it is known that interesting transformations of entanglement occur j22|. Results 
not shown here indicate circumstances under which such intra-dimer entanglement can be 
broken if a dimerized state is superposed with a non-dimerized, completely random state. In 
particular this results in dimer density matrices being very close to Werner states [23] and 
therefore results in the entanglement of dimers vanishing when the random state component 
is more than 2/3. In general the effect of superposition on entanglement has been studied 
vigorously j9[, as well as entanglement in the RVB states has been explored (Hf. The results 
discussed in this paper may add in some small measure to the understanding of entanglement 
in such contexts. 

Real states are often obtained as eigenstates of time-reversal symmetric systems and the 
discussion here will be of relevance to such systems, for instance to spin chains that have 
time-reversal. Clearly when the states are complex the above approach to finding optimality 
conditions does not work. In the case when k = 2 the optimality condition is found to be 
the same as that for real states, namely < riir2%/r\ 2 < 1 (except now r\2 is not restricted 
to the reals). Thus it is conceivable that as this is the unit interval in the complex plane, 
the measure of optimal states is zero. However of course, there could be an infinity of these, 
for instance all real states are possible candidates. 



IV. DISCUSSIONS AND SUMMARY 
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Appendix: Evaluation of an integral for the optimal fraction f 2 

To recall, the integral is 

r 4 

h = / 6^22)9(^2 -r u r 22 )P({x i })P({yJ) \{d Xi dy h 
J i=i 



(A.l) 



where = (ipi\<J y <E> Cy\4>j) , = 1,2 are two real and normalized 4- vectors, and the 
measures P are the uniform measures in Eq. ( |T5l) . First, the particular Pauli matrix o y 
that appears can be replaced by the other Pauli matrices. In particular the a z matrix being 
diagonal, offers a simpler look. This replacement is quite easily seen to be equivalent to 
some 45° rotations of the original variables. 

Also dropping the constraint on the product rnr 2 2 will be useful. If the resultant integral 
is denoted as /, then it is shown below that f 2 is simply / — 1/2. Next, it is proven that as 
far as / is concerned, the two 4-vectors can be taken to be orthogonal. Decompose say \ip 2 ) 
along the vector j^i) and one orthogonal to it: 

|V*)=cos(0)|Vi)+sm(0)b7>> ( A -2) 

where (77 1 ^1 ) = 0. No additional phases are involved as the states are all real. A straight- 
forward calculation shows that 

r\ 2 - r u r 22 = sin 2 (6>) ((Vik* ® ^zVlf ~ (vWz ® o^K^iK ® <Tz\*l>i)) • (A. 3) 

The quantity within the brackets is precisely the same combination as in the L.H.S., except 
that instead of the vectors are orthogonal. As sin 2 (#) has a constant positive sign, this 
proves that we can consider the pairs, to begin with, as being orthogonal. In other words the 
sign of the combination r\ 2 — r u r 22 is invariant under the Gram-Schmidt orthogonalization 
process. 

The additional constraint of the vectors being orthonormal introduces an additional Dirac 
delta function term in the measure. Writing the fraction / as a ratio f n / fd, the numerator 
/„ and the denominator fd are given by 

/» = ^ / 6 ( r ?2 - rnr 22 ) 5 (j^ x iVi J 6 A ~ 1 J * Vi ~ 1 J dx dy, (A.4) 

and 

f*= J / 5 fe J 5 (E " x ) 5 ^ " x ) dx ^ ( A - 5 ) 

All the sums are from 1 to 4 and dx dy is the Euclidean eight dimensional volume ele- 
ment. The factor 8/tt 2 is introduced for later convenience alone. To be further explicit the 
combination r\ 2 — r u r 22 = 

(xxyi - x 2 y 2 - x 3 y 3 + x 4 y 4 ) 2 - (x\ -x 2 2 -xj + x\){yl - y\ - y\ + y\). (A.6) 



Introducing a series of transformation to various two-dimensional polar coordinates (in the 
(£1,24) pair, the (22,23) pair, etc., as well as in the resulting radii) and performing two 



12 



delta function integrals corresponding to the normalizations results is: 

f n = 2 6 [(cos(a) cos(/3) cos(0) — sin (a) sin(/3) cos(0)) 2 — cos(2a) cos(2/?)] 

J (A. 7) 

S [cos (a) cos(/3) cos(#) + sin (a) sin(/3) cos(0)] sin(2a) sin(2(3)dad(3d9d(j), 

and a corresponding integral for fd, only without the Heaviside theta function constraint. 
Here a, /3 G [0, 7r/2] while 0, <fi G [0, 27r]. As a check, the integral can be easily done without 
either the Heaviside theta or the Dirac delta functions to give 8tt 2 , which is exactly the 
factor that follows from the normalization of the two Dirac delta normalization constraints; 
see Eq. (|T5|) . if one takes into account the factor 8/tc 2 that is introduced in Eq. (1A.4j) . 

Introducing variables v and u as cos(a) cos(/3) cos(#) ± sin(a) sin(/3) cos(0) respectively, 
allows the delta function integration over the v variable to be performed and results in: 

, 6 {u 2 — cos(2a) cos(2/3)) du da d(3 

Jn= I — i I (A.8) 

u 2 I u 2 




4cos 2 (a) cos 2 (/3) V 4 sin 2 (a) sin 2 (f3) 

Notice that the given range of a and 0, [0,7r/2] 2 can be divided into four equal squares, 
such that the G function constraint is effective only in [0,7r/4] 2 and in [7r/4, 7r/2] 2 , as 
cos(2a) cos(2/9) is negative elsewhere. The range of the u integration is restricted depending 
on a, (3. Taking the range [0, vr/4], the contribution from it is denoted 11 for lower-lower, 
we have: 

/•vr/4 fir/4 pui j 

fS= da d/3 n du , ==, (A.9) 



where 



J/3 Ju I U 2 U 2 




1 - 



4c 2 c 2 V 4s 2 s 2 



Po = sin 1 a/1/2 - s 2 , M = ^c 2a c 2 i3, and u x = 2s a sp, 



and s a stands for sin(a) etc.. The limits of the integration are such that the function 
constraint is satisfied as well as the square-roots are real numbers. The denominator fraction 
can be similarly split up, and in fact f l J differs from the above in that /?o = as well as 
u = 0. It is also not hard to see that f% u (uu for upper-upper) is same as and similarly for 
fd. In the other two regions, as the constraint is not operational, and because of symmetry, 
it follows that: ff = ff = fn = fn = f ul - It then follows that 

An evaluation of the three corresponding integrals is carried out numerically and results in 
fl = 0.116850..., f ul = 1/2, f l J = 0.285391.... It is remarked that standard softwares could 
not evaluate the integrals symbolically, however the numerical results are sufficient to give 
the fraction / = 0.785398..., which is tt/4 to 1 part in 10 6 . It is also easy to similarly see 
that f l j = (n — 2)/4 and from consistency/" = {n 2 — 8)/16. 

Returning to the integral in Eq. (I A . 1 1) for the fraction f 2 we see that there are two O 
constraints, while we have considered only one, namely the second one. If the first constraint 
0(^11^22) alone is present, it is easy to see that the integral is 1/2. This also follows from 



13 



the fact that r n and r 22 are both independent and uniformly distributed. We also state 
parenthetically without proof that r 12 is distributed according to the semi-circle distribution. 
Now, if r n r 2 2 is negative, then surely r\ 2 — rnr 2 2 is positive, which is 50% of the time. Thus 
the fraction of cases when r n r 22 is positive and r\ 2 — r n r 22 is positive is / — 1/2, which 
is precisely the required fraction / 2 . Thus we have evaluated / 2 and presented sufficient 
evidence that it is actually (ir — 2)/4. It is not clear if it is only a coincidence that this is 
also precisely The evaluation presented here may, by far, not be the "optimal" one, but 
is the best the authors could come up with. 
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